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ABSTRACT 

We report observations of two consecutive transits of the warm super-Earth 
exoplanet GJ 1214b at 3.6 and 4.5 /im with the Infrared Array Camera instrument 
on-board the Spitzer Space Telescope. The two transit light curves allow for the 
determination of the transit parameters for this system. We find these paremeters 
to be consistent with the previously determined values and no evidence for tran- 
sit timing variations. The main investigation consists of measuring the transit 
depths in each bandpass to constrain the planet's transmission spectrum. Fix- 
ing the system scale and impact parameters, we measure Rp/R-k = 0.1176to;oQ09 
and O.llGSloiooos 3.6 and 4.5 /im, respectively. Combining these data with 
the previously reported MEarth Observatory measurements in the red optical 
allows us to rule-out a cloud- free, solar composition (i.e., hydrogen-dominated) 
atmosphere at 4.5 a confidence. This independently confirms a recent finding 
that was based on a measurement of the planet's transmission spectrum using 
the VLT. The Spitzer, MEarth, and VLT observations together yield a remark- 
ably flat transmission spectrum over the large wavelength domain spanned by 
the data. Consequently, cloud-free atmospheric models require more than 30% 
metals (assumed to be in the form of H2O) by volume to be consistent with all 
the observations. 

Subject headings: planetary systems — methods: observational — techniques: 
photometric 
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Introduction 



GJ 1214b is a transiting; super-E arth exoplanet orbiting a nearby M4.5V star ev- 
ery 1.6 days (ICharbonneau et al.l 120091) that was discovered by the MEarth Observatory 



( iNutzman fc Charbonneaij|2008l : llrwin et al.ll2009l ). The first data for this planet suggested 
that it has a mass of 6.55 ± 0.98 and a radius of 2.68 ± 0.13 R(q, and subsequent follow- 
up transit hght curve observations a r id analysis have confirm ed the initial radius estimate 
(IBerta et al.ll2010l : ICarter et al.ll2010l : iKundurthy et al.ll2010t ). Due to degeneracies between 
possible interior and atmospheric ma sses and compositions, a diverse range of models can re- 
produce G J 1214b's mass and radius (IRogers fc Seageiil2010bl : iNettelmann et al.ll2010l). This 



problem had been previously identified for planet s in the super-Earth regime (jSeager et al. 



20071 : lAdams et al.l l2008l : iRogers fc Seagerl l2010af ). Nevertheless, GJ 1214b is larger than 



expected for its mass assuming a purely solid composition and mature age; this implies that 
GJ 1214b most probably possesses a significant gas envelope. An important question then 
is, what are the characteristics of this atmosphere? 

The transiting nature of GJ 1214b fortunately prov ides the opportunity to study its at- 
mosphere in transmission (e.g. ICharbonneau et al.ll2002l ). In particular, the overall composi- 
tion of the planet's atmosphere (e.g., hydrogen-dominated or metal-rich) can be investigated 
using transm ission spectroscopy obse r vations through the determiri ation of the atmospheric 
scale height (IMiller-Ricci et al.l 120091 : iMiller-Ricci fc Fortneyl l2010l ) . Because the depth of 
strong absorption features in the transmission spectrum is proportional to the atmospheric 
scale height, such observations constrain the mean molecular weight of the atmosphere, and 
thus its bulk composition. Obtaining constraints on the compositi on of GJ 1214b's atmo- 
sphere can give clues to the planet's origins and interior structure (IRogers fc Seagerl l2010bl : 
Nettelmann et allboiol ). 



In this Letter, we present transit observations of GJ 1214b obt ained with Spitzer that 
yield measurements of the planet's transmission spectrum. Recently, iBean et al.l ( l2010l ) also 
presented a transmission spectrum for the planet in the red optical from VLT observations. 
We use our Spitzer data along with the VLT and MEarth data to determine the composition 
of GJ 1214b's atmosphere, and we discuss the results in the context of the possible interior 
structure models and origins for the planet. 



2. Warm-Spitzer Observations 



GJ 1214b was observe d for approxirnately 4.5 h during each of two consecutive transits 
with Warm-Spitzer /IKAC (IWerner et al.ll2004j : lFazio et al.ll2004l ) at 3.6 and 4.5 /xm (program 
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ID 542). The observations occurred on 26 and 27 April 2010, with the first transit observed 
at 3.6 yum. Observations in each visit were gathered in sub-array mode (32 x 32 pixels) with 
an exposure time of 1.9 s per frame, which resulted in a 2.0 s cadence. These yielded 109 
images per visit, with 64 frames per image and a total of 6976 frames. Nearly one-fourth of 
these exposures were taken during the 53 min planetary transit. 

The method we used to produce photometric time se ries in each channel from the Basic 



Calibrated Data (BCD) is described in lDesert et al.l (120111 ). It consists of finding the centroid 



position of the stellar point spread function (PSF) and performing aperture photometry. 
The aperture size was optimized to minimize the final residuals in each channel, although we 
noticed that the S/N does not vary significantly with the aperture radii. The final aperture 
sizes were set to 5.5 and 2.5 pixels at 3.6 and 4.5 /xm, respectively. The background level 
for each frame was determined by fitting a Gaussian to the central region of a histogram 
of counts from the full array. We find that the background varies within the sets of 64 
frames, but that it remains globally constant between images and for the whole time-series. 
The 1.9 s exposures gave a typical signal-to- noise ratio (S/N) of 340 and 280 per individual 
frame at 3.6 and 4.5 /xm, respectively. The noise is 20% and 10% larger than photon noise 
at 3.6 and 4.5 /im, respectively. We used a sliding median filter to select and trim outliers 
greater than 5 a, which correspond to 1% of the data. We also discarded the first half- 
hour of observations, which are affected by a significant telescope jitter before stabilization. 
The final number of photometric measurements used are 6194 and 5998 for the light curves 
gathered at 3.6 and 4.5 /im, respectively. The two final raw time series are presented in 
Figure [TJ 



3. Determination of the transit parameters 
3.1. Analysis of the Warm-Spitzer light curves 



We used a transit light curve model multiplied by instrumental decorrelation functions 
to measure the t ransit parameters and their uncertainties from the Spitzer data as described 
in iDesert et al.l (1201 1[). We computed th e transit light curves with the IDL transit routine 



OCCULTQUAD from iMandel &: Agoll ( 120021 ). which depends on five parameters: the planet-to- 
star radius ratio Rp/Ri,, the orbital semi-major axis to stellar radius ratio (system scale) 
a/R^,, the impact parameter b, the time of mid transit Tc, and a linear limb-darkening co- 
efficient Ci. In the case of late M-dwarfs, such as GJ 1214, the limb-darkening is uncertain 
because it strongly depends upon stellar physical parameters that are not known accurately. 
Consequently, we opted to determine the limb- darkening coefficients for our light curve mod- 
eling from the data themselves. We assumed that the limb-darkening is well-approximated 
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by a linear law ( lClaretlll998l ) at infrared wavelengths. 



The Spitzer /YRKC photometry is known to be systematically affected by the so-called 
pixel-phase effect (see e.g., Charbonneau et al. 2005). This effect is seen as oscillations 
in the measured fluxes with a period of approximately 70 min (period of the telescope 
pointing jitter) and an amplitude of approximately 2% peak-to-peak. We found the centroid 
position of the PSF with GCNTRD, from the IDL Astronomy Library 0. We decorrelated our 
signal in each channel using a linear function of time for the baseline (two parameters) and a 
quadratic funct ion of the PSF posit ion (four parameters) to correct the data for each channel 
as described in iDesert et al.l ( 120111 ) . 



We performed a simultaneous Levenberg-Marquardt least-squares fit ( lMarkwardtll2009[ ) 
to the data to determine the transit and instrumental parameters (11 in total). The errors on 
each photometric point were assumed to be identical, and were set to the rms of the residuals 
of the initial best-fit o btained. To obtain an estimate of the correlated and systematic 
errors (jPont et al.ll2006l ) in our measure ments, we use the r esidual permutation bootstrap, or 
"Prayer Bead", method as described in lDesert et al.l ( l2009l ). In this method, the residuals of 
the initial fit are shifted systematically and sequentially by one frame, and then added to the 
transit light curve model before fitting again. We allow asymmetric error bars spanning 34% 
of the points above and below the median of the distributions to derive the 1 a uncertainties 
for each parameters. 



3.2. Results 

We performed two separate analysis of the Spitzer data, and the results are given in 
Table [1] For the first analysis, we fit the data without any constraints and allowed all the 
transit parameters to be unique for each of the two Spitzer channels. Our values for all 
the parameters are in agreeme nt with the previously reported resul t s from extensive follow- 



up observations of GJ 1214b flBerta et al.ll2010l : ICarter et al.ll2010l : iKundurthv et al.l 12010 



Sada et al.l l201o[) and th e adop ted parameters used to determine the optical transmission 
spectrum by iBean et al.l ( 20101 ). Furthermore, the mid-transit times from the Spitzer ob- 
servations are precisely measured wit h an accurac y of ± 7 s and do not deviate significantly 
from the revised linear ephemeris of iBerta et al.l (120101 ). This strongly suggests that the 
transiting parameters we derive are robust and well constrained at this point. 

The agreement between all the available data sets justifies a second analysis of the Spitzer 



^ http : //idlastro . gsf c . nasa . gov/homepage . html 
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data, where we fixed the b and a/i?^ values to those used by iBean et aLl (120 lOf ) to measure 
the optical transmission spectrum while still determining unique Rp/ i?^ for each channel and 
collapsing along the dominant degeneracies bet ween these paramet ers. We also reanalyzed 
the MEarth data for five transits presented by iBerta et al.l (120101 ) with these parameters 
fixed. This enables us to combine the Spitzer photometry with the VLT spectrum and the 
MEarth photometry to constrain GJ 1214b's transmission spectrum. 

Adopting identical b and a/R-^, the MEarth data yield Rp/R^, = O.llSltaoooy, and the 
average radius ratio value in the VLT spectrum is Rp/ Ri, = 0.1165 with a typical uncertainty 
of 0.0006. For the Spitzer data, we find Rp/R^ = 0. 11761^-^^^g and 0.1 1631^-^^;^^ at 3.6 and 



4.5 fim, respectively (see Tabled]). As described by iBean et al.l ( 120101 ) . the VLT spectrum 
is featureless, and we find that the transit depths derived from the Spitzer, MEarth and 
VLT measurements are also all in agreement within 2 a level. The linear limb-darkening 



(Claret 


2000; 


Sine 


2010) 



The period of the telescope jitter leading to the pixel-phase effect has the same timescale 
of the transit duration (1 hour). Therefore, stronger degeneracies exist between the transit 
depths, the linear limb-darkening coefficients, and the corrections of the intra-pi xel sensitivity 



than is observed in transit data for planets with longer transit durations (e.g. iDesert et al. 
20091 ). This leads to somewhat larger error bars than had been expected, although fixing the 
limb-darkening coefficients only reduces the the error bars on the transit depths by 20%. 

Stellar variability can be responsible for time and wavelength dependent transit depth 
variations ( iDesert et al.ll201ll ). Since the sets of observations we aim to compare have been 
secured at different epoc hs, it is critical to e nsure that GJ 1214's variability does not affect the 
measured radius ratios. iBerta et al.l ( l2010l ) analyzed three years of photometric monitoring 
in various bandpasses and found the stellar rotation period of GJ 1214 to be 52.7 ± 5.3 
days, and the photometric variability in the red optical to be 1%. This suggests that the 
transit depths vary by less than 0.01% and 0.004% (in absolute) in the MEarth and Spitzer 
bandpasses, respectively, which is below the amplitude errors of the radius ratios. This is also 
further supported by multiple transit observations of this object in various bandpasses cited 
above. Consequently, the comparison between the transit depth derived from observations 
secured at different epoch can be done with no corrections. 
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4. Discussion 

4.1. The nature of GJ1214b's atmosphere 

As already mentioned in Sect. I3.2[ the wavelength dependent measurements of Rp/Ri, 
secured by Spitzer, MEarth, and VLT are all in agreement within the 2 a level. In principle, 
the lack of features in the wavelength dependent planet radius can be interpreted as no 
evidence for the presence of a planetary atmosphere seen in transmission. This means that 
either this planet has no atmosphere or its depth variations are below our level of sensitivity. 
However, our current knowledge of the planetary mass and radius rules out the absence 
of an atmospheric envelope since it is required b y interior structure models of GJ 1214b 



(IRogers fc Seageii l2010bl : iNettelmann et al.l 120101 ). Therefore, we assume that GJ 1214b 



possesses a gas envelope in the rest of this paper. 

Figure [2] shows the Spitzer, MEarth, and VLT Rp/ R^ measurements compared to cloud 
free tran smission spectroscopy rnodels for GJ 1214b with various gas compositions calcu- 



lated by iMiller-Ricci fc Fortneyl (120 lOl ). The atmospheric models solve for the wavelength- 
dependent deposition of flux, as a function of depth, and make no assumption regarding 
albedo. We assume heat redistribution over the full planet, which results in an equilibrium 
temperature of 555 K. None of these models include cloud formation since abundant conden- 
sates are not expected to form at the temperatures expected for the planet's atmosphere. We 
consider two main types of atmospheric models. The first type includes hydrogen-dominated 
atmospheres with metallicities ranging from solar to 50 times solar, and also a solar com- 
position model with methane artificially removed. The later is considered because methane 
molecules are susceptible to photolysis since they have large photo-dissociation cross sections 
in the UV. Furthermore, methane does not re-form very quickly through chemical reactions 
at the temperatures that are expected for the upper atmosphere of GJ 1214b. The second 
type of models considered are mixtures of molecular hydrogen and water (in the form of 
vapor) ranging from 10 to 100% water by volume. All the spectral models were computed at 
high resolution and were fiux weighted integrated over the bandpasses of the measurements 
to enable calculation of the goodness of fit. We apply a scaling factor to the models to 
give the best fit to the measurements. 

As a first test, we compare the models to only the MEarth and the Spitzer data. We 
find that these data on their own exclude hydrogen-dominated atmosphere models with 
solar metallicity at 4.5 a confidence level. This is due to a strong broad band methane 
absorption feature at 3.6 /xm which disa ppears in the methane-free model as seen in Figure [2J 



This is consistent with the findings of iBean et al.l (l2010l ) using only the VLT spectrum 



The hydrogen-dominated atmospheres would have a large scale height and lead to large 
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spectral variations that are not observed in the data. The difference between the radius 
ratios measured in the two Spitzer bandpasses, ARp/R^ = 0.0013^q;qq}2 limits the range of 
possible atmospheres because of the methane feature. Of the cloud-free hydrogen-dominated 
atmospheres, only those without methane are consistent with the combined Spitzer and 
MEarth data. 

The combination of MEarth^ VLT, and Spitzer observations provides stronger con- 
straints on the models, in particular thanks to the spectral information from the VLT data. 
The combined data set rules out the model with solar metallicity at 7 a confidence level and 
models with up to 50 times enhanced metallicity at 6.3 a. This is due to the sensitivity of the 
da ta from the VLT b andpass to spectral features of water, which has been noted previously 
by iBean et al.l (|2010[ ). Including the VLT data also enable us to rule out the methane- free 
model at 5 cr confidence. 

The flat transmission spectrum of the combined data set suggests that the scale height 
must be small for cloud-free atmospheres to produce spectral features with an amplitude 
lower than 0.02%. An atmosphere dominated by chemical elements heavier than hydrogen 
and helium is therefore required. Although we cannot discriminate between heavy elements, 
water in its vapor state is expected to be a do minant species in various model scenarios 
( iRogers fc Seagerll2010bl : iNettelmann et al.ll2010l ). The MEarth, VLT, and Spitzer observa- 
tions together are consistent within 3 a with atmospheric models containing more than 10% 
of water by volume. Models with more than 30% of water are consistent within 1 a. 

An alternative explanation for the featureless transmission spectrum of GJ 1214b indi- 
cated by the MEarth, VLT, and Spitzer data could be the presence of high altitude clouds or 
hazes. A high altitude cloud deck would reduce the amplitudes of the expected spectral fea- 
tures because the stellar light would be transmitted through a smaller atmospheric layer than 
in the case of a cloud-free at mosphere. The proposed temperature-pressure profiles of the 
Miller-Ricci fc Fortneyl (|2010[ ) models do not cross the condensa tion curves of known equi- 
librium species that are expected to be abundant (lFortneyll2005l ). Furthermore, clouds and 
hazes particles have wavelength dependent opacities which seem at odd with the remarkably 
fiat transmission spectrum observed. However, more theoretical work is needed since, for 
example, a distribution of particle sizes or photochemical processes, although unconstrained 
at this point, could modify the transmission signature of this atmosphere. Hydrocarbon 
haze formed through photochemical processes and methane depletion by similar processes 
remain possibilities to explain the broad band observations. Consequently, we only consider 
the simple case of a cloud and haze-free atmosphere at chemical equilibrium in the following 
section. 
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4.2. The atmosphere as a diagnostic of the planet's origins and interior 



Rogers &: Seagerl (l2010bl ) propose three major distinct possible origins for the gas layer 
concomitant with interior structure models for GJ 1214b: accretion from the primordial 
protoplanetary nebulae, sublimation of ices, or outgassing of rocky material. Our observa- 
tions exclude the primordial atmosphere scenario (with possible exceptions such as a hazy 
methane-free atmosphere) because this would yield a hydrogen-dominated composition. At- 
mospheres arising from sublimation or outgassing are consistent with the data. 

If sublimation of ices is the atmospheric origin scenario, it would suggest that the planet 
has a high concentration of water in its interior. This implies that the planet formed beyond 
the snow line of its host star's protoplanetary disk, and that it migrated inward to its 
current orbital position. After the planet migrated closer to its star, the ice forming the 
bulk sublimated and led to a metal-dominated atmosphere as it is currently observed. There 
are two possible scenarios to explain the hydrogen-poor envelope in this case. Either the 
planet never attained a mass sufficient to accrete or retain large amounts of primordial gas 



(H-rich) from the protoplanetary disk (jRafikovll200(:), or the mass fractions o f light elements 



have been diminished through atmospheric escape. iRogers fc Seaged (l2010bl ) point out that 
an atmosphere with a small scale height, such as suggested by our observations, is more 
robust against atmospheric escape than a hydrogen-rich envelop, and that the cumulative 
loss of mass from a water-dominated atmosphere would not significantly affect the overall 
character of the planet. Interestingly, iNettelmann et al.l (120101 ) show that a water-rich planet 
that completely lacks hydrogen in the atmosphere requires an implausibly large water-to- 
rock ratio. Instead, H/He/H20 envelopes with water mass fractions between approximately 
50 and 85% are favored. This picture is in agreement with the MEarth, VLT, and Spitzer 
observations, which suggest a water mass fraction in the atmosphere of more than 50%. 

The alternative scenario for the origin of GJ 1214b's atmosphere would require a gi- 



lope fe.e. Pollack et al. 


1996), or outgassing during formation and evolution or a period of 


tectonic activity 


(Schaefer & Feelev 


2007; 


Elkins-Tanton & Seaeer 2008: Schaefer & Feelev 


2009: Kite et al. 


2009h. 


Rogers & Seaeei 


( 


2010bf) discussed this scenario in the context of an 



atmosphere composed of purely hydrogen on the grounds that heavy elements like water or 
carbon dioxide can not be out-gassed in sufficient quantity to form an atmosphere that can 
account for GJ 1214b's large size on their own. However, they did not consider the possibility 
of a hybrid atmosphere with significant amounts of both hydrogen an d heavy elements. While 



heavy elements alone can not be out-gassed in sufficient quantities, lElkins-Tanton fc Seagei 



(120081 ) have shown that signific ant amounts of water can be out-gassed from chondritic plan 



etesimals. Given this, and that INettelmann et al.l (120101 ) have indicated a mixed H/He/H20 
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envelope can successfully be used to model the planet, a hybrid composition atmosphere 
from outgassing seems to be a plausible scenario. 
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Fig. 1. — Spitzer transit light-curves observed in the two currently available IRAC band- 
passes at 3.6 (top) and 4.5 /xm (bottom). Top panels : raw transit light-curves. The red 
solid line corresponds to the best fit models which include the time and position instrumental 
decorrelations as well as the model for the planetary transit (see details in Sect. [2]). Middle 
panels : corrected and normalized transit light-curves with their best fits. The data are 
binned in 100 seconds intervals (50 points). Bottom panels : residuals. 
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Table 1. GJ 1214b parameters and uncertainties from Spitzer observations. 



Wavelength 


Rp/ -R* 


b 




Cl 


Tc (BJDsTo) " 




O-C (See.) b 


3.6 fim 
4.5 fim 


1176+00010 

U.il/D_o 0010 
„ ,,-,,+0.0011 
U.ii<i_o 0015 


10+01^ 

n 2';+o-ii 

'-'•^"^-0.16 


15 46+0-13 
1= 2t;+0.39 

10.ZO_o 5g 


r, r,,o + 0.054 

U.^i»_0.050 
n 146+0 088 
U.i4D_o 051 


55312.633877tgg^I 
55314.21423ltgg^I 


05 
05 
05 
05 


-13+12 

_n+i2 

11-12 


3.6 fj,m 
4.5 fim 


n 1176+0 0008 

U.ii/D_o 0009 

n 1163+0 0010 

U.ilD^>_o 0008 


fixed to 0.27729 
fixed to 0.27729 


fixed to 14.9749 
fixed to 14.9749 


„ 2'3q+0.044 
u.z,3y_0 

n sn5+o io8 

U.OUO_o 067 


2455312.633085t7:g^ 
2455314.213440lg;3^ 


-05 
-05 
-05 
-05 


_in+i2 

11-12 

-9+12 
^-12 



''Times are given as Baryeentric Julian Dates in the Barycentric Dynamical Time system iEastman et al.ll2010h . 

''Observed minus Calculated (O-C) using the ephemeris from l|Berta et al.ll20ld ) with Period (days)=l. 58040490 ± 0.00000033 
and To (BJDsdt)=2454966.525042 ± 0.000065. 



- 11 - 




Fig. 2. — Comparison of the measured wavelenKth dependent planet- to-star radius ratios to 
transmission spectroscopic models from lMiller-Ricci fc Fortneyl ( 120101 ) . The two radius ratios 
obtained from the Spitzer observations are the black filled circles with their 1 cr error bars. 
The results from MEarth and VLT observations are shown as black square and triangles, 
respectively. The continuous lines correspond to the best -scaled and smoothed transrn ission 
spectra expected from model of GJ 1214b atmosphere's (IMiller-Ricci &: Fortneyll2010l ). The 
red continuous line corresponds to a transmission spectra with solar composition (H/He 
rich), the red dotted lines is the same spectra but with metallicity enhanced by a factor of 
50. The green line is a model with solar composition but with no methane and the blue line 
is a model with 100% water vapor in the atmosphere. The open circles are the flux weighted 
integrated models in the Spitzer bandpasses. The dotted black lines at the bottom of the 
plot correspond to the instrumental bandpasses. The hydrogen rich model (red continous 
line) is ruled out a 7 a level by the combined set of observations. 
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